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Abstract. Measurements of NO x and ozone performed during the NOXAR project are 
compared with results from the coupled chemistry-climate models 
ECHAM4.L39(DLR)/CHEM and GISS-model. The measurements are based on flights 
between Europe and the East coast of America and between Europe and the Far East in the 
latitude range 40°N to 65°N. The comparison concentrates on tropopause altitudes and 
reveals strong longitudinal variations of seasonal mean NO, of 200 pptv. Either model 
reproduced strong variations 3 km below but not at the tropopause, indicating a strong 
missing NO x or NO } sink over remote areas, e.g. NO x to HN0 3 conversion by OH from 
additional OH sources or HN0 3 wash-out. Vertical profiles show maximum NO x values 2- 
3 km below' the tropopause with a strong seasonal cycle. ECHAM4.L39(DLR)/CHEM 
reproduces a maximum, although located at the tropopause with a less pronounced 
seasonal cycle, whereas the GISS model reproduces the seasonal cycle but not the profile’s 
shape due to its coarser vertical resolution. A comparison of NO x frequency distributions 
reveals that both models are capable of reproducing the observed variability, except that 
ECHAM4.L39(DLR)/CHEM shows no very high NO x mixing ratios. Ozone mean values, 
vertical profiles and frequency distributions are much better reproduced in either model, 
indicating that the NO x frequency distribution, namely the most frequent NO x mixing ratio, 
is more important for the tropospheric photochemical ozone production than its mean 
value. 

Both models show that among all sources, NO x from lightning contributes most to the 
seasonal cycle of NO x at tropopause altitudes. The impact of lightning in the upper 
troposphere on NO x does not vary strongly with altitude, whereas the impact of surface 
emissions decreases with altitude. However, the models show significant differences in 
lightning induced NO x concentrations, especially in winter, which may be related to the 
different treatment of the lower stratospheric coupling between dynamics and chemistry. 
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1. INTRODUCTION 


Global numerical simulations have been widely used in the past years to evaluate the influence 
of anthropogenic emissions on climate. Tropospheric gas phase chemistry modeling focussed 
strongly on ozone because of its role as a greenhouse gas. It has a pronounced sensitivity in 
radiative forcing peaking at the tropopause altitude (Wang and Sze, 1980). In addition, model 
studies suggest stronger climate sensitivity than well-mixed greenhouse gases (Ponater et al., 
1999). Recent results suggest that global tropospheric ozone increase since pre-industrial times 
lead to a radiative forcing of 0.3 to 0.6 W m' 2 (WMO, 1999) enhancing significantly the global 
mean forcing of well-mixed greenhouse gases of 2.25 to 2.5 W m' 2 (IPCC, 1996). Furthermore, 
numerical simulations were used to assess the inf Aence of the increasing air traffic which, leads 
to significant NO,-emissions at flight altitude (IPCC, 1999). 

The abundance of reactive nitrogen compounds, such as NO, N0 2 , and N0 3 (=NO,), has changed 
significantly since pre-industrial time (Lelieveld and Thompson, 1999). Major contributors to the 
concentration of upper troposphere NO* are lightning, aviation and NO, emissions at the surface, 
which largely consist of energy consumption, soils, vegetation and biomass burning. Ozone 
production, and therefore climatic impacts, depend on the NO, concentration (Fishman and 
Crutzen, 1979; Fishman et al., 1979). Coupled chemistry-climate models, used to assess climate 
impacts from NO, emissions, need to be able to simulate realistically NO, and ozone. 

However, the level of confidence of the magnitude of the increase of tropospheric ozone is low 
(IPCC, 1996), because the knowledge is mainly based on numerical simulations. The physical 
and chemical processes, which determine global tropospheric ozone concentrations, are very 
complex. They need to be simplified in global numerical models mainly because of limitations of 
modern computer facilities. Comparisons of the NO,-fields at 200 hPa calculated by ten state of 
the art models show pronounced disagreement with respect to the absolute concentrations and 
the horizontal variations (Brasseur et al., 1998). The first systematic comparison between the 
available NO, measurements with numerical simulations was made by Emmons et al. (1997). 
However, NO, measurements at the tropopause altitude have been rather limited in the past 
because they were restricted to aircraft campaigns. 

The NOXAR (Nitrogen OXides and ozone measurements along Air Routes) measurements now 
enable an extensive model validation with regard to upper troposphere NO, and ozone. These 
measurements were performed from a passenger aircraft during a one-year period ending in May 
1996. The permanently installed measurement system recorded chemical and meteorological 
data on 270 return flights between Europe and destinations in the United States of America and 
Far East, and is therefore one of the largest data sets for upper troposphere NO, (Brunner, 1998). 

The present investigation has two primary aims: First, an extensive comparison of NOXAR 
measurements with the climate-chemistry models ECHAM4.L39(DLR)/CHEM and the GISS 
model will give a better understanding of upper troposphere dynamical and chemical processes. 
Second, sensitivity studies with both models provide an interpretation of the NOXAR 
measurements in terms of the origin of upper troposphere NO, and ozone. 

The ECHAM4.L39(DLR)/CHEM and the GISS model have similar concepts for the description 
of the meteorology which is described by the physical part of the model. It provides detailed 



information for the chemistry module about the state of the atmosphere, like cloud cover and 
structure, precipitation, temperature, and radiation, necessary for the calculation of photolysis 
and reaction rates, NO, sinks (e.g. wet deposition) and sources (e.g. lightning). This is different 
to many chemical transport models, which use prescribed meteorological fields as provided by 
e.g. ECMWF and climatological values for cloud processes. The tropospheric chemical schemes 
are quite similar in either model. Differences between the models, which are important in the 
context of this study, include vertical resolution, different lightning parameterizations and the 
treatment of the lower stratosphere. ECHAM4.L39(DLR)/CHEM includes a description of the 
lower stratosphere dynamics and chemistry, which is not explicitly treated in the GISS model but 
based on prescribed fixed fields (see next section). In this comparison very similar NO, emission 
strengths are used. 


2. DATA-SET AND MODELS 


2.1. NOXAR measurements 

The NOXAR measurement system was permanently installed on a Swissair B-747 airliner. The 
instrumentation package consisted of commercial analyzers for the measurements of NO, 
(N0 2 +N0 3 ) (chemiluminescence) and ozone (UV-Photometer). Calibrations were automatically 
initiated in regular time steps. Detection limits of 23, 45, and 500 pptv were achieved for 2- 
minute mean concentrations of NO, (N0 2 +N0 3 ), and ozone, respectively. The accuracies are 8, 
1 1 and 6%, respectively. (N0 2 +N0 3 ) concentrations were calculated from daytime NO and 0 3 
measurements assuming a photostationary steady state. Measured values were only used for 
nighttime flights when available. Due to instrumental problems, no (N0 2 +N0 3 ) was measured 
between May and November 1995. Average daytime NO, concentrations based on measured and 
on calculated (N0 2 +N0 3 ) agreed within 11%. Preference was given to the calculated values 
because of the rather low precision of the (N0 2 +N0 3 ) measurement. Details are given in Brunner 
(1998). Data, which were identified to be influenced by young aircraft plumes were excluded. 

Meteorological parameters were recorded from the aircraft data system. The NOXAR 
measurements cover mid-latitudinal flight tracks, mainly between 40°N and 65°N and from 
Eastern parts of North America to Central Europe and to China (Figure 1). Dias-Lalcaca et al. 
(1998) and Brunner (1998) give detailed descriptions and include also more southern flights to 
Bombay and Hong Kong. On those flights, the airplanes were usually flying several kilometers 
below the tropopause so that they are not included here. 

Figure 1: Area of measurements 

2.2. ECHAM4.L39(DLR)/CHEM model 

The dynamical part of the coupled climate-chemistry model ECHAM4.L39(DLR)/CHEM is the 
ECHAM4.L39(DLR) model. It represents an advancement of the operationally used general 
circulation model ECHAM4 (Roeckner et al., 1996) regarding the vertical resolution, which has 
been increased from 19 to 39 levels from the surface to the top layer, centered at 10 hPa (30 km) 
(Land, 1999; Land et al., 1999). This model version has been designed especially for upper 
troposphere and lower stratosphere processes. It includes state of the art parameterizations of 
radiation, cloud, surface and soil processes, gravity wave drag, vertical diffusion, convection, 
and semi-Lagrangian transport for water vapor, cloud water and chemical species. The model 
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generates its own climate, which depends only on prescribed sea surface temperatures (Gates, 
1992) and the atmospheric C0 2 concentration for 1990 and either interactively calculated (online 
chemistry coupling) or prescribed (offline chemistry coupling) concentrations of the radiatively 
active gases H 2 0, 0 3 , CH 4 , N 2 0, and CFCs. A detailed validation of the ECHAM4.L39(DLR) 
model is given in Land (1999). 

The chemistry module CHEM (Steil et al., 1998) has been coupled to ECHAM4.L39(DLR) 
(Hein et al., 2000). It includes stratospheric heterogeneous and homogeneous ozone chemistry 
and tropospheric N0 x -CH 4 -C0-0H-0 3 -chemistry with 107 photochemical reactions and 37 
species. Surface emissions of NO x due to industry, soils, and biomass burning are based on 
Benkovitz et al. (1996), Yienger and Levy (1995) and Hao et al. (1990), respectively. Aircraft 
NO x emissions are based on the DLR2 data set (Schmitt and Brunner, 1998), NO x emissions by 
lightning are calculated interactively according to Price and Rind (1994) (see Table 1). 
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Table 1: Brief comparison of the coupled chemistry-climate models 

ECHAM4.L39(DLR)/CHEM and GISS. 

2.3. GISS model 

The GISS H' ('two prime') GCM (Hansen et al., 1997a) was developed from the GISS GCM 
version II which is described fully by Hanson et al. (1983). GISS II’ is a gridpoint model having 
a 4x5 degrees horizontal resolution and 9 sigma-pressure hybrid levels in the vertical. 
Climatological, monthly sea surface temperature data based on the Atmospheric Modeling 
Intercomparison Project (AMIP) 1979-1993 and ocean ice coverage are incorporated. The 
radiation scheme determines radiative heating in the short- and longwave spectral regions for all 
important radiatively active species in the troposphere and stratosphere (Hansen et al., 1983). 
The II' version employs several improvements over model II, including new convection (Del 
Genio and Yao, 1992), prognostic cloud (Del Genio et al. 1996), planetary boundary layer 
(Hartke and Rind, 1997) and ground hydrology (Rosenzweig and Abramopoulos, 1997) 
treatments, as well as improvements to the albedo calculation over land (Rossow and Zhang, 
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1995), ocean (Gordon and Wang, 1994) and sea ice (Hansen et al., 1997b). Tracers are advected 
using the quadratic upstream method (Rind and Lemer, 1996), which produces little numerical 
noise and realistically maintains strong tracer gradients. A detailed description can be found in 
Rind and Lemer (1996). 

The GISS model includes 10 chemical tracers (O,, NO,, N 2 O s , HN0 3 , H 2 0 2 , CH 3 OOH, HCHO, 
H0 2 N0 2 , CO and H 2 0) using the 'family' approach, where appropriate with 52 reactions for 24 
species (DeMore et al., 1997). Peroxyacetyl nitrate (PAN) and explicit non-methane hydrocarbon 
(NMHC) chemistry is omitted. The 'FASTJ' photolysis scheme, after validation with results from 
the Photochemical Model Intercomparison (PhotoComp, Olson et al., 1997), is employed. N 2 0 5 
hydrolysis on sulfate aerosol is calculated according to Dentener and Crutzen (1993) using 
monthly mean modeled sulfate surface area taken from Koch et al. (1999). Wet and dry 
depositions are adapted from Wang et al. (1998). The lightning parameterization generates NO, 
interactively with the convection scheme (Price et al., 1997). Chemical emissions are prescribed 
as monthly-mean values and include CO (biomass, equivalent isoprene and industrial) and NO, 
(industrial, biomass burning, soils, and aircraft). The GISS climate-chemistry model is 
documented extensively and validated against observational data in Shindell et al. (2000). 

2.4. Experimental set-up 

For the comparison of the model results with the NOXAR measurements (sections 3 to 5) two 
numerical simulations performed with the ECHAM4.L39(DLR)/CHEM and one performed with 
the GISS model were taken into account (Table 2). The ECHAM4.L39(DLR)/CHEM model was 
run in two modes: on-line and off-line. In the on-line mode feedbacks from the chemistry module 
to radiation via changes in ozone, water vapor, methane, nitrous dioxide, and CFCs are allowed, 
whereas in the off-line mode those species are prescribed for the radiation module with 
climatological values. Two and one year spin-up times are employed for the 
ECHAM4.L39(DLR)/CHEM and the GISS model, respectively. Since the coupling of the upper 
troposphere and lower stratosphere is included in the ECHAM4.L39(DLR)/CHEM model, 
whereas the GISS model has there its upper boundary, and since the chemical lifetime of NO, is 
larger in the lower stratosphere than in the upper troposphere, a longer spin-up time has to be 
taken into account for ECHAM4.L39(DLR)/CHEM. After the spin-up time inter-annual 
variability occurs, but no clear trend is visible. 

In addition to the control runs shown in Table 2, perturbation simulations were performed with 
both models in order to estimate the impact of various NO, sources on the NO, and ozone 
distribution in the upper troposphere measured during the NOXAR project (section 6). For each 
NO, source (aircraft, lightning, surface) simulations have been performed with both models, 
which neglect each source in turn. For the GISS model the impact of the upper boundary has also 
been investigated by setting the NO, upper boundary equal to zero. The 
ECHAM4.L39(DLR)/CHEM model was run for these simulations in the off-line mode and the 
GISS model in the on-line mode. Five and two years simulations were performed including a 
spin-up time of two and one years for the ECHAM4.L39(DLR)/CHEM and the GISS model, 
respectively. 
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Table 2: Experimental set-up for the control simulations for the 

ECHAM4.L39(DLR)/CHEM and the GISS models. 

2.5. Model differences 

Regarding the dynamical part of the GISS and the ECHAM4.L39(DLR)/CHEM model, main 
differences occur in the vertical resolution of 9 and 39 layers and the tracer transport schemes, 
quadratic upstream and semi-Lagrangian method, respectively. Numerical diffusion caused by 
the coarser resolution of the GISS model may be compensated by the transport scheme, which 
maintains strong tracer gradients. 

The chemical schemes and NO, sources and sinks of either model are quite similar. However, 
the GISS model includes an additional NO, sink (N 2 O s hydrolysis on sulfate aerosol) and takes 
cloud to cloud and cloud to ground flashes into account and applies a C-shape NO, emission 
profile, whereas ECHAM4.L39(DLR)/CHEM applies a uniform mixing for cloud to ground 
flashes. Although the total amount of NO, released into the atmosphere is similar in both models, 
the GISS model releases significantly more in the upper troposphere due to the differences in the 
parameterization . 

An important difference is also the coupling of dynamical-chemical processes in the upper 
troposphere and lower stratosphere. The GISS model has an upper boundary for NO, and 0 3 at 
1 10 hPa, which excludes transport of NO, from the tropical upper troposphere to the lower mid- 
latitude stratosphere. The ECHAM4.L39(DLR)/CHEM includes stratospheric chemistry and 
accounts for those processes. However, this difference may only be of minor importance for the 
control experiments, since those processes are implicitly included in the boundary condition, but 
they definitely affect the sensitivity tests. 

The GISS model includes a variety of online diagnostics, like budget calculations, wave analysis, 
etc. This enables a more detailed analysis compared to the ECHAM4.L39(DLR)/CHEM model. 

2.6. Methodology: Tropopause determination 

The tropopause defies a sharp transition from the troposphere to the stratosphere, where ozone 
strongly increases. To avoid systematic differences in the comparison between model and 
observational data, the comparison of NOXAR measurements with model data was performed 
for the tropopause region. This region represents air masses at the tropopause ±25 hPa for the 
measurements and the ECHAM4.L39(DLR)/CHEM model data and the 7 th level of the GISS 


6 




i 


model (200 hPa), where the GISS model has its temperature minimum regarding mean vertical 
profiles. The height of the tropopause and its pressure were calculated for every measurement 
and model time step. Data from the European Centre for Medium-Range Weather Forecasts 
(ECMWF) were used to calculate height and pressure of the dynamic defined tropopause with a 
constant value of the potential vorticity (2 PV-Units) and to calculate the distance from the 
measurements to the tropopause. The height and pressure of the thermal defined tropopause were 
calculated for ECHAM4.L39(DLR)/CHEM model data. The use of two different approaches to 
estimate the height of the tropopause does not limit the results since the thermal defined 
tropopause can be well represented by a constant value of potential vorticity in the region of the 
measurements (Grewe and Dameris, 1996). 


3. LONGITUDINAL VARIATIONS 

Figure 2 shows the variability of NO, mixing ratios along the flight tracks between the eastern 
part of the United States of America and China. In all seasons, large values are observed over 
America including a seasonal cycle with 150 pptv in winter and 350 pptv in summer. Minimum 
values were measured over the eastern part of the North Atlantic, varying between 100 pptv and 
150 pptv in winter and summer, respectively. Over Eurasia enhanced values were seen over 
Central Europe (10°E to 15°E) in every season and also over parts of Asia, especially in summer, 
where NO, exceeds 350 pptv. This means that a large longitudinal variability is observed with an 
amplitude of 1 00 pptv in winter and 200 pptv in summer. 

Figure 2: Longitudinal variations of NO, 

Both models, the ECHAM4.L39(DLR)/CHEM (on- and off-line version) and the GISS model, 
show some longitudinal variation (Figure 2, dotted and solid thin lines), but do not reproduce the 
observed magnitude. ECHAM4.L39(DLR)/CHEM simulates enhanced NO, volume mixing 
ratios over America and Europe with an amplitude of 30 pptv in winter and 60 pptv in summer. 
But it misses totally the enhanced Siberian summertime values (see below). Generally, mean 
values are too low especially over the continents by up to 200 pptv, thus the land-ocean contrast 
is less pronounced. In contrast, the GISS model has a zonal uniform NO, distribution at mid 
latitudes in all seasons with a strong seasonal cycle. The inter-annual variability is much larger in 
the GISS than in the ECHAM4.L39(DLR)/CHEM model. 

The discrepancies, revealed by the comparison, can be restraint to two possible causes: 

(1) The upper tropospheric NO , lifetime is too long, i.e. NO , sinks are too weak or even 
missing, especially over the North Atlantic, where a minimum in the observed NO, mixing ratios 
occurs in all seasons. The inclusion of PAN-chemistry in the models would decrease even further 
the land-sea contrast in NO, as it would be converted into PAN near the sources, i.e. over the 
continents, and produced over remote areas, i.e. the sea. Instead, the inclusion of an additional 
OH source, like acetone, which has similar source regions as NO,, but a longer lifetime, would 
provide a source of OH in remote areas (Arnold et al., 1997). This leads to a NO, sink over the 
North Atlantic due to the reaction N0 2 + OH + M -4 HN0 3 + M. 
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(2) N0 X production from HN0 3 photolysis may be too strong, which implies either to large 
mixing or again a to small total sink of NO y at tropopause altitude. The tropopause region is a 
sharp transition between the tropospheric circulation and the stratospheric circulation, which is 
dominated by zonal structures. At 300 hPa, i.e. approximately 3 km below the tropopause, both 
models show strong latitudinal variations in NO x of 35 to 150 pptv and 200 to 400 pptv in the 
ECHAM4.L39(DLR)/CHEM and GISS model, respectively. The horizontal structures of NO y 
are similar to those of NO x at each altitude for either model. Both models show strong changes of 
the simulated patterns of NO x and NO y from 300 to 200 hPa, indicating changes of dynamical or 
chemical processes, in contrast to observations. At 200 hPa, both models overestimate the mean 
zonal wind by 5 to 10 ms' 1 and show to low transient wave energy, which accounts for horizontal 
patterns, in mid-latitudes at tropopause levels (Land, 1999). This leads to a more intense mixing 
of NO x and NO y at tropopause altitudes. However, the deviations from observations are not large 
enough to levei off NO x and NO y variations. Therefore, this leads again to a missing sink. A sink 
of NO y can only be wash-out of HN0 3 , suggesting that modeled wash-out processes me 
underestimated in both models at mid-latitude tropopause region. For mid-latitudes observed 
vertical profiles of precipitation and cloud water are not available. First results of observations in 
the tropics (Tropical Rainfall Measuring Mission, TRMM) suggest that ocean storms produce 
precipitation water at higher altitudes compared to land storms (Del Genio, personal 
communication, 2000). Since the GISS model simulates precipitation water at lower altitudes 
and without any land-sea contrast, one could speculate that this also occurs at mid-latitudes. This 
would lead to an underestimation of the washout of HN0 3 especially over the North Atlantic. 

A comparison of the zonal mean chemical NO x loss rate at 200 and 300 hPa in the latitude band 
40°N to 60°N in July reveals for the GISS model large differences of 1.4x10V and 22x10V, 
respectively (see also Figure 4). This leads to a change of the process, which controls NO x in the 
GISS model: chemistry at 300 hPa and dynamics at 200 hPa. This supports reason (1) that the 
NO x sink is underestimated at 200 hPa. However, it does not exclude (2). 

The upper troposphere source for NO x in the ECHAM4.L39(DLR)/CHEM seems to be 
underestimated, whereas in the GISS model this source is overestimated in spring, summer, and 
fall. Thus ECHAM4.L39(DLR)/CHEM shows mean NO x values, which are typically North 
Atlantic, i.e. around 100 pptv with a small seasonal cycle, whereas the GISS model reflects 
typically land situations, such as over the eastern part of the U.S.A. or Europe, with higher 
values and a strong seasonal cycle. 

Although ECHAM4.L39(DLR)/CHEM shows enhanced NO x values over America and Europe in 
summer, the model fails to simulate enhanced values over Siberia. Measurements of lightning 
events performed with the Optical Transient Detector (OTD, Christian et al., 1996) give evidence 
for strong lightning in that area in summer, which locally has the same order of magnitude as the 
maximum values observed in tropical areas. Since the model simulates much less lightning there 
(not shown), a significant NO x source is missing in that area, which is probably the main reason 
for that discrepancy. 


4. VERTICAL DISTRIBUTION 

Figure 3 shows the mean vertical NO x profiles for each season. Observed maximum values are 
located 2-3 km below the tropopause. A clear seasonal cycle is observed with values 
approximately twice as high in spring and summer than in fall and winter. 
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Figure 3: Vertical Profile of NO, 


As already stated above, the ECHAM4.L39(DLR)/CHEM model shows too low NO, values for 
all seasons at the tropopause altitude. However, it is able to simulate a local maximum at the 
tropopause, i.e. 2 to 3 km higher than observed. The measured seasonal cycle with low values 
(75 pptv) m winter and high values (150 pptv) in summer is quite well reproduced regarding the 
summer to winter ratio: it shows twice as high values in summer than in winter. Due to the 
coarse vertical resolution, the GISS model only shows a local maximum at the tropopause in 
July, but simulates a reasonable seasonal cycle. 

The more detailed diagnostics of the GISS model enables a closer look at the discrepancy of the 
ocation of the summertime NO, maximum. At 200 hPa and in the latitude range 40°N to 60°N 
and 175 hPa to 260 hPa, lightning is the dominant NO, source (10.5 kg s ') and meridional and 
vertical transport the dominant NO, sink (9.7 kg s 1 ) (Figure 4). Below, at 260 hPa to 400 hPa, 
convection (14 kg s' ) is a more important source for NO, than lightning (5.8 kg s 1 ) and transport 
(2.6 kg s ). Chemistry is the only important sink (18 kg s '). Although the total NO source is 
smaller at 200 hPa than at 300 hPa the NO, mixing ratios at 200 hPa exceed those at 300 hPa 
because of the larger NO, residence time at 200 hPa, which in term depend on the sink terms. In 
section 3 we suggested that a strong NO, sink is missing at tropopause levels, e.g. 200 hPa. This 
would imply a decrease of the location of NO, maximum, supporting the findings in section 3. 

Figure 4: NO, Budgets 

Ozone profiles are viewed in Figure Figure 5. Measurements show troposphere values of 40 to 
50 ppbv in winter and 60 to 70 ppbv in summer and a sharp increase at the tropopause for all 
seasons. Although the models obviously have problems in representing the longitudinal NO 
variations, the vertical ozone distribution agrees very well with the observations. Even the sharp 
gradient at the tropopause is well reproduced by the ECHAM4.L39(DLR)/CHEM model in 
winter and spring and is only slightly smeared in other seasons. Despite its coarser vertical 
resolution, the GISS model also simulates a sharp gradient at the tropopause. 

Figure 5: O, Profiles 


5. TEMPORAL VARIABILITY 

The frequency distribution of NO, at tropopause levels, derived with the NOXAR measurements 
(Figure 6) shows broad distributions with peak values between 65 and 85 pptv. NO x mixing 
ratios of less than 50 pptv occur clearly less often in summer than in other seasons And instead 
high NO, mixing ratios of more than 200 pptv occur more often in summer.’ 
ECHAM4.L39(DLR)/CHEM (on-line version) simulates the most frequent NO, mixing ratios 
between 40 and 60 pptv and the GISS model in the range of 40 pptv to 1 10 pptv. Thus the most 
frequent NO, mixing ratios are slightly too low and too confined’ in the 

ECHAM4.L39(DLR)/CHEM model and too low in winter and too high in summer in the GISS 
model. 


The comparison indicates that the NO, variability is quite well simulated by both models with 
two exceptions. (1) The main differences between measurements and 
ECHAM4.L39(DLR)/CHEM data are the underestimation of high (more than 125 pptv) NO, 
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mixing ratios. (2) The main differences between measurements and the GISS model is its too 
strong seasonal cycle with too low values in winter and too high values in summer. 

Thus the comparison of the frequency distributions explains quite well the differences in the 
mean values shown in Figure 2. Mean modeled NO, mixing ratios are too low in 
ECHAM4.L39(DLR)/CHEM, not only because of a systematic shift in the frequency 
distribution, but mainly due to a lack of the occurrence of high values of more than 100 pptv in 
winter and 150 pptv in summer. This means that large-scale NO, plumes with concentration 
increases of typically 300 pptv and up to 5 ppbv as described in Brunner et al. (1998) are not 
observed in the model. However, it does not explain the land-sea contrast viewed in Figure 2, 
since the measured frequency distributions at North Atlantic sites still show occurrences of air 
masses with high NO, content, although less frequent than over land (not shown). The GISS 
model shows a shift of the most frequent occurrence of NO, mixing ratios to higher values in 
summer (solid line). This leads to an overestimate of the NO, shown in Figure 2. 

The ozone production depends non-linear on the NO, concentration. Therefore a realistic 
simulation of the frequency distribution, namely the most frequent NO, mixing ratio, is more 
important for a correct modeling of the total ozone production than mean NO, values. This is 
indicated by the model results: Most frequent NO, mixing ratios (Figure 6) are simulated fairly 
well and also ozone mixing ratios (Figure 5), but not the mean NO, values (Figure 2). 

Figure 6: Density function NO„/0 3 

The density function of measured ozone at tropopause levels (Figure 6) shows the most frequent 
occurrences at around 50 ppbv in winter and autumn and at higher values (75 ppbv) in summer 
and spring. Also the variability is significantly increased in spring and summer. Figure 6 shows a 
systematic shift of the frequency distribution to higher ozone values in the 
ECHAM4.L39(DLR)/CHEM model. This may not be over-interpreted, since this is related to the 
location of the tropopause, which is a sub-grid parameterization. As shown in Figure 5, ozone 
has a sharp gradient at the tropopause in the measurements as well as in the modeled data. A 
vertical shift of the tropopause by some hundred meters will result in dramatic changes of the 
ozone values. Comparison of the width of the density functions yields a good agreement with 
lower variability in autumn and winter and higher in spring and summer for the 
ECHAM4.L39(DLR)/CHEM model. The GISS model shows much broader distributions and 
only small variations between the seasons. That means that the GISS model shows a large 
seasonal cycle for NO, but not for ozone and the ECHAM4.L39(DLR)/CHEM model for ozone 
but not NO,. 


6. ORIGIN OF NO x AND OZONE IN THE UPPER 

TROPOSPHERE 

To determine the origin of upper troposphere NO, and ozone, several sensitivity simulations 
were performed for both models. In each simulation one NO, source (lightning, aircraft or 
surface) was excluded. For the GISS model an additional simulation with an upper boundary 
condition for NO, (equal to zero) was done. The simulation lengths are 2 and 5 years including 1 
and 2 years spin-up for the GISS and the ECHAM4.L39(DLR)/CHEM models, respectively. 
Results are summarized in Table 3 and Table 4 for the mid-latitude Northern Hemisphere upper 
troposphere and they are compared to other model studies. Lightning (5% to 35%) and surface 
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NO, sources (9% to 30%) contribute with a similar amount to the upper tropospheric NO, in 
winter and are doubled in summer. The contribution due to aircraft emissions shows a large 
uncertainty among the studies, 11% to 60% in winter and half these values in summer. The 
stratospheric impact is less then 30% in winter and less than 20% in summer for all models 
except the GISS model. 

This shows that in winter the low vertical mass exchange by large-scale ascent and convection 
diminish the importance of surface as well as of lightning NO, emissions and increases the 
importance of aircraft emissions. In the GISS-model, the importance of the stratosphere is quite 
large, 60% in winter and 34% in summer. Figure 4 shows that the downward transport of 
stratospheric NO, is the dominant source for the upper troposphere NO, in the GISS model in 
January, whereas it is only a minor contributor in July. Other model studies (Krauss et al., 1996; 
Kohler et al., 1997; Brasseur et al., 1998) s ’owed a much reduced stratospheric impact of 2% to 
30% and 5% to 25% in January and July, respectively, which suggests that results from the GISS 
model are influenced by the upper boundary. Consequently, the impact of other NO x sources is 
underestimated, so that the GISS model lies at the lower side of the estimated contributions of 
the surface and aircraft NO x source for January. 

The large uncertainties of the importance of lightning for the upper troposphere NO x burden, 
revealed in the comparison, can be attributed to the parameterization of the lightning NO x source. 
Although the total amount of lightning produced NO x is similar in each model simulation (5 and 
4.7 TgN a' 1 ), the vertical distribution is not. In the GISS model vertical profiles derived from 
measurements (Pickering et al., 1998) are applied, leading to much higher emissions in the upper 
troposphere than below. Other model studies distributed the NO, emitted per lightning flash 
uniformly with height. The effect is clearly seen, as the GISS model shows a much higher impact 
of lightning in the upper troposphere in summer than other model studies. Applying C-shape 
lightning emission profiles (Pickering et al. 1998) to the ECHAM4.L39(DLR)/CHEM model 
shows a 70% to 75% contribution of the lightning NO, source to the upper troposphere mid 
latitude NO, burden, confirming these findings. 


January 

E4C 

GISS 

Kr96 

Ko97 

Br98 

Br98* 

Total 

Lightning 

18 

18 

<10 

5 

30-35 


5-35 

Surface 

25 

9 

30 

25 

10 


9-30 

Aircraft 

51 

11 

60 

57 

15 

13-26 

11-60 

Stratosphere 

- 

60 

<10 

2 

15-30 


2-60 


Table 3: Contribution [in % relative to the control simulation] of various NO, sources to the 
atmospheric NO, content at tropopause levels and in the area of the NOXAR measurements (see 
Figure 1) for the ECHAM4.L39(DLR)/CHEM and GISS model. Results of the model studies 
Krauss et al. (1996) (Kr96), Kohler et al. (1997) (Ko97) and Brasseur et al. (1998) (Br98) are 
added for the northern hemisphere mid latitudes (30°N to 60°N) at 200 hPa. Brasseur et al. 
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(1998) concentrated on the effect of aircraft emissions and presented results from an inter- 
comparison of 4 models (Br98*). 


July 

E4C 

GISS 

Kr96 

Ko97 

Br98 

Br98* 

Total 

Lightning 

36 

56 

20 

16 

40-50 


16-56 

Surface 

31 

26 

40 

56 

20 


20-56 

Aircraft 

17 

5 

10 

16 

25 

31-32 

5-32 

Stratosphere 

- 

34 

20 

5 

10 


5-34 


Table 4: As Table 3 but for July. 

Figure 2 showed an increase of measured NO, in the upper troposphere over land in summer 
compared to winter. To attribute this increase to distinct sources we assume that the strength of 
the lightning, aircraft, and surface NO, emissions and the stratosphere flux account for the 
changes in upper troposphere NO,. The relative importance r’’ of each NO, source P (lightning, 
aircraft, surface, and stratosphere) is calculated by r 1 ’ = [(C S -F)-(C W -P W )] / { X Jp [(C s -P s ) - (C w - 

P w )l } , where C and P denote the NO, mixing ratio at northern upper troposphere mid-latitudes in 
the control and perturbation simulations, respectively. The indices s and w refer to summer and 
winter values. At mid-latitude tropopause altitudes the NO, to NO y ratio decreases in the 
ECHAM4.L39(DLR)/CHEM model from 25% to 15% from 30°N to 60°N and is around 30% 
for the GISS model for both, summer and winter season. So that the ratio remains mainly 
unchanged between summer and winter. The SONEX and POLINAT2 campaigns showed 
similar ratios for fall (Singh et al., 1999; Talbot at al., 1999; Ziereis et al., 2000). At mid latitudes 
and tropopause levels, chemical conversion affects winter and summer values in the same 
manner and therefore can be neglected for this approach. 

Table 5 summarizes the results. All models show that at least 70% of the increase in summer 
compared to winter can be attributed to lightning or surface NO, emissions. However, some 
models, including the ECHAM4.L39(DLR)/CHEM and the GISS-model, identify lightning, 
while others identify surface emissions as the main source for the summer to winter difference. 
In the ECHAM4.L39(DLR)/CHEM and the GISS-model, the NO, contributions due to lightning 
and surface emissions show the same horizontal pattern, reflecting that NO, concentrations at 
tropopause altitudes from both sources are related to convection (not shown). Jeker et al. (2000) 
analyzed observational data derived during the POLINAT2 and SONEX campaign. They found 
high NO, concentrations at the tropopause and could attribute them to lightning events by using 
backward trajectories and additional satellite data. However, in some cases surface NO, still 
could contribute by a large portion. 
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E4C 

GISS 

Kr96 

K697 

Br98 

Total 

Lightning 

67 

55 

27 

26 

48 

26-67 

Surface 

46 

26 

49 

83 

25 

25-83 

Aircraft 

-23 

2 

-3 

-17 

29 

-17-29 

Stratosphere 

GO) 

17 

27 

8 

-3 

-3-27 


Table 5: Contribution [%] of different NO, sources to the increase of upper troposphere NO, in 
northern mid latitudes from winter to summer. For the ECHAM4.L39(DLR)/CHEM model the 
same transport characteristic regarding stratospheric NO, as for the ECHAM3 (Kohler et al., 
1997) is assumed. 

Figure 7 shows the contributions of the different NO, sources to the NO, and ozone mixing ratio 
as a function of height. The GISS model shows an aircraft contribution to the atmospheric NO, 
content of 10% to 20% in the region from the tropopause to four kilometers below in winter and 
summer. The contribution due to lightning increases between those heights from around 15% in 
winter to 40%-50% in summer, which leads in summer to the local maximum at the tropopause 
(Figure 3). The importance of surface NO, emissions decreases with height from 50% to 10% in 
winter and 45% to 30% in summer. 

Figure 7: Changes in the vertical profile of NO, and ozone due to distinct NO, sources 

The ECHAM4.L39(DLR)/CHEM model shows a similar behavior regarding the importance of 
the surface NO, sources in the upper troposphere. Maximum impact of surface sources occur 1 to 
2 km below the tropopause in summer (30% to 40%) and 6 km below the tropopause in winter 
(80%). Since ECHAM4.L39(DLR)/CHEM also includes stratospheric chemistry and dynamics, a 
propagation of upper troposphere NO, into the lower stratosphere via the Hadley cell circulation 
is enabled. This leads to an accumulation of lightning and aircraft NO, emissions in the lower 
stratosphere, where the NO, lifetime is much longer than in the troposphere due to washout of 
HNOj. The stratospheric NO, can return to the troposphere and increase NO, and ozone there, 
Maximum impacts of lightning and aircraft NO, emissions are calculated to be in the lower 
stratosphere and at tropopause regions, respectively. 

The similar shape and magnitude of the impact of surface NO, sources also shows that the two 
different parameterizations for convective mass transport used in either model lead to similar 
results. Another possibility to evaluate vertical mixing away from the surface are Radon 222 
( 222 Rn) simulations, since 222 Rn is mainly emitted over land and has a short half-life of 3.8 days. 
The comparison of 222 Rn simulations performed with both models (Land, 1999; Rind and Lemer, 
1996) show very similar global distributions of 2 2 2 Rn, confirming that 
ECHAM4.L39(DLR)/CHEM and the GISS model simulate vertical mixing in a similar way. It 
indicates that the model differences regarding NO, mainly arise from the upper boundary 
condition and lightning parameterization. 
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Figure 3 showed a strong summertime increase of measured NO, going from the lower 
stratosphere and upper troposphere to approximately 4 km below the tropopause. Both models 
suggest that this increase predominantly results from lightning and surface NO, emissions. Their 
importance varies with altitude. The importance of the lightning emissions increases, whereas the 
importance of the surface emissions decreases with altitude. At tropopause altitudes it is mainly 
lightning which causes the NO, increase in summer, whereas 4 km below the tropopause surface 
emissions are the dominant factor. 

Figure 7 (right panel) shows the ozone changes due to various NO, sources calculated with the 
ECHAM4.L39(DLR)/CHEM (top) and the GISS model (bottom). The aircraft induced ozone 
changes in the lower stratosphere are relatively small compared to the high ozone background. 
Both models show a several percent change of ozone due to aircraft emissions in the upper 
troposphere. Surface NO, emissions have a si nilar effect on ozone in both models. Their impact 
is largest in the middle troposphere and decreases with decreasing altitude, i.e. toward the 
tropopause. However, ECHAM4.L39(DLR)/CHEM shows a larger impact in winter, whereas the 
GISS model implies this in summer. Since the ECHAM4.L39(DLR)/CHEM model accounts for 
stratospheric chemistry and since it shows a strong increase of NO, in the lower stratosphere due 
to lightning, therefore the impact of lightning on ozone (15% to 20%) is much stronger than 
shown by the GISS model (5% to 10%). 

Figure 8 shows the impact of each NO, sources on the frequency distributions of NO, and ozone. 
The ECHAM4.L39(DLR)/CHEM model shows that aircraft emissions not only increase the NO, 
burden in the upper troposphere but also its temporal variability. Aircraft emissions are even the 
main contributor to the variability in winter. The frequency distribution broadens, and especially 
NO, mixing ratios between 40 and 100 pptv occur oftener. The effect of lightning and surface 
emissions is similar, but not as pronounced. In summer, all sources contribute significantly to the 
NO, variability. 

Figure 8: Changes in the frequency distribution of NO, and ozone due to distinct NO, 
sources 

The GISS model shows reasonable NO, frequency distributions. In winter lightning, surface and 
aircraft NO, emissions change the occurrence of low NO, values (less than 200 pptv) only little. 
Instead, stratosphere-troposphere exchange, namely the upper boundary condition, controls the 
occurrence of high NO, values of more than 100 to 150 pptv. In summer, the impact of the upper 
boundary is of much less importance. Instead, lightning contributes significantly to the higher 
summertime NO, values, as well as to its variability. 

Both models agree that lightning contributes most to the NO, variability among all NO,-sources. 
However, surface emissions are also of importance in both models and contribute significantly to 
the variability. The models show very different behavior concerning aircraft emissions. This is 
related to the treatment of stratospheric transport processes (see discussion above). 

The ozone distribution at tropopause levels is strongly determined by stratospheric ozone. The 
impact of tropospheric NO, sources is much smaller (see discussion above). Consequently the 
ozone variability is also hardly affected. Only lightning may change the frequency distribution 
significantly, as shown by the ECHAM4.L39(DLR) model. In winter, lightning leads to a 
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broadening of the ozone frequency distribution, whereas in summer only a shift to higher values 
is modeled. 


7. CONCLUSIONS AND DISCUSSION 

Aircraft measurements of NO x and ozone derived within the Swiss NOXAR project have been 
extensively compared with two chemistry-climate models. The comparison has been 
concentrated on the tropopause region on flight routes between Europe and the eastern part of 
America and between Europe and the Far East. For the first time chemistry models have been 
intensively validated with regard to upper troposphere mean NO x mixing ratios and its 
variability. 

The observations show large longitudinal variations of NO x . Although the models show regional 
patterns of NO x at tropopause altitudes, they do not simulate the strength of the observed 
longitudinal variations. This clearly indicates a missing NO x or NO, sink, especially over the 
North Atlantic. The ECHAM4.L39(DLR)/CHEM model shows too low NO x mixing ratios with 
only small longitudinal variations, suggesting that the NO x , which is released in the upper 
troposphere is underestimated, and it will be even more underestimated when taking into account 
an additional NO x sink. However, sensitivity tests showed a doubling of NO x mixing ratios, when 
applying C-shape instead of uniform NO x emission profiles for flashes, although neither the 
horizontal distribution nor total amount of the NO x lightning emissions was changed. 

At 300 hPa, i.e. approximately 3 km below the tropopause, both models are able to simulate 
strong longitudinal variations indicating that the observed differences are confined to the 
tropopause altitude. 

We examined the reasons for the observed to modeled NO x differences and found two 
possibilities. Either an additional NO x sink is missing or HN0 3 washout is underestimated by 
either model. For a more detailed examination, budgets of HO x , NO x and NO y are required, 
including information about cloud water and precipitation, which are not available. The models 
do not include acetone, which would provide an additional OH source over remote areas, a 
potential candidate for a missing NO x sink (Arnold et al., 1997). However, Brune et al. (1999) 
showed that the main sources of HO x in the northern mid-latitude upper troposphere are water 
vapor, formaldehyde, and peroxides, with acetone being a smaller contributor. On the other side, 
there are some indications that HN0 3 washout could be underestimated: First results of the 
Tropical Rainfall Measuring Mission showed larger cloud water amounts due to storms over 
oceans than over land. The GISS model does not show these differences and simulates cloud 
water amounts shifted to lower altitudes (Del Genio, personal communication, 2000). If these 
findings were also be confirmed for the mid-latitudes, it would give reasons that HN0 3 washout 
is underestimated with a distinct land-sea contrast. 

The online diagnostics of the GISS model showed distinct differences between chemical- 
dynamical interactions at 200 hPa and 300 hPa. At 200 hPa, i.e. tropopause altitudes, NO x is 
mainly controlled by transport, since the chemical net destruction rate is relatively small. The 
main sinks are meridional and vertical transport, and lightning is the dominant source. At 
300 hPa, i.e. 3 km below, NO x is controlled by chemistry, with convection as the main source. It 
indicates again that the transition from the troposphere to the stratosphere is not correctly 
simulated and that the total NO x sink is underestimated at tropopause altitude. 
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Although effort has been taken to avoid systematic differences between measurements and 
model data regarding the sampling method (e.g. comparison relative to tropopause altitude), 
there are still some possibilities for systematic sampling differences. Model data are taken 
uniformly over time and space. On the other hand, the aircraft measurements were sampled over 
various times during day and night, but not uniformly. The aircraft entered the Shannon control 
area between 6 and 8 UTC and 13 and 17 UTC for east- and westbound flights (Brunner et al., 
2000). However, since day and nighttime measurements are included, systematic differences 
should be small. A more important difference in the sampling of observational and model data 
may be the fact that the flight routes avoid jet stream areas on west-bound flights, and storm 
systems entirely. A classification of the NO* data with regard to the weather situation and the 
flight tracks would be necessary. However, the large longitudinal variations are quite unlikely to 
result solely from sampling differences. 

Mean vertical profiles of NO* show an E-shape, with maximum values at the surface, below the 
tropopause, and in the stratosphere, and with minimum values in the free troposphere and 
lowermost stratosphere. This shape could be reproduced by ECHAM4.L39(DLR)/CHEM model 
and partially also by the GISS model. However, observed maximum values tend to be some 
kilometers below the tropopause, whereas the models simulate it directly at the tropopause. 
Again this is probably related to the simulation of the troposphere-stratosphere transition. It has 
been shown for the GISS model that the residence time of NO* at 200 hPa exceeds those at 
300 hPa due to a dramatic change in the chemical sink. It indicates again a missing or 
underestimated NO, sink. Radon 222 experiments yield very similar results for either model 
(Land, 1999; Rind and Lemer, 1996). Since Radon 222 has a short have-life, this indicates that 
transport, namely convection and large-scale ascent is simulated similar in either model. It points 
out that both models show a shift of the NO* maximum to higher altitudes for the same reasons. 

However, simulated ozone distributions agree well with observations. Even observed sharp 
vertical gradients at the tropopause were realistically reproduced. This indicates that the 
horizontal distribution of NO* and its mean value is of minor importance for the ozone 
distribution at tropopause levels. Instead the frequency distribution and the most frequent NO* 
mixing ratio are of larger importance, since the local net ozone production depends on these 
values rather than on the mean values. These NO* frequency distributions are reasonable 
simulated by both models with two exceptions. ECHAM4.L39(DLR)/CHEM is missing high 
NO* mixing ratios (greater than 1 50 pptv in winter and greater than 200 pptv in summer) and the 
GISS model has too many occurrences of high NO* mixing ratios in summer. Both models are 
able to reproduces reasonably the most frequent NO* mixing ratios. 

The ECHAM4.L39(DLR)/CHEM model was run in to modes, with and without interactively 
coupling of the chemistry with dynamics. Both model versions showed similar results, justifying 
the use of the off-line version for sensitivity tests. 

The ECHAM4.L39(DLR)/CHEM model (off-line version) and the GISS model were used to 
estimate the impact of individual NO* sources on the NO* and ozone distribution. They showed 
that the seasonal changes in NO* are mainly due to lightning NO* emissions. The impact of 
lightning emissions on the NO* and ozone distribution in the upper troposphere is mostly 
constant or slightly increasing with altitude, whereas the impact of surface emissions strongly 
decreases. Lightning also contributes most to the NO, variability in summer. In winter, aircraft 
emissions contribute most to the NO* variability in the ECHAM4.L39(DLR)/CHEM model, 
whereas stratosphere to troposphere exchange is the dominant contributor in the GISS model. 
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Although there are several qualitative agreements between the models, also discrepancies occur: 
Due to the different treatment of the chemistry and dynamics in the lowermost stratosphere, the 
effect of aircraft and lightning emissions on the NO x distribution and the effect of lightning NO x 
emissions on the ozone distribution differs significantly. ECHAM4.L39(DLR)/CHEM shows a 
much larger impact of aircraft emissions on the NO x distribution at the tropopause especially in 
winter and a much larger effect on ozone due to lightning. 

The horizontal patterns of the upper troposphere NO x distribution clearly mirrors the emissions, 
i.e. maximums over industrial regions and regions with high convective activity. This indicates 
that the ratio of the chemical NO x lifetime to the dynamical lifetime is small enough to allow for 
regional patterns. However, the comparison to observations shows clearly that the chemical NO x 
lifetime is overestimated at tropopause altitudes. The estimated importance of individual NO x 
sources should not be affected thereby, since emphasis is givm to relative changes. 
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9. Figures 

Figure I: Area covered by NOXAR measurements. Thin lines indicate the land-sea mask in T30 
resolution, as used by the ECHAM4.L39(DLR)/CHEM model. 

Figure 2: Seasonal NO x volume mixing ratios (pptv) in the tropopause region (± 25 hPa relative 
to the local tropopause) and as a mean of all latitudes covered by the NOXAR measurements. 
The very heavy solid, heavy solid, the solid and the dotted lines represent NOXAR 
measurements, the online version, the offline version of ECHAM4.L39(DLR)/CHEM model, 
and the GISS-model, respectively. The inter-annual standard variation is added for the 
models. The number of NOXAR measurements per grid box (i.e. per 3.75°) is added. 

Figure 3: Mean vertical NO, profiles (pptv) relative to the tropopause for each season. The very 
heavy solid, heavy solid, the solid and the dotted lines represent NOXAR measurements, the 
online version, the offline veision of ECHAM4.L39(DLR)/CHEM model, and the GISS- 
model, respectively. The inter-annual standard variation is added for the models. The number 
of NOXAR measurements per grid box (i.e. per km) is added. The high spring at 2 to 4 km 
below the tropopause are cause by individual large plumes, and may not represent the 
climatological mean value. 

Figure 4: NO x Budgets for the GISS model for the 7 th (200 hPa) and 6 th (300 hPa) level in the 
latitude range 40°N to 60°N for January (italic) and July (bold). NO x exchange due to 
horizontal and vertical transport, convection, chemical processes, lightning and aircraft 
emissions are given in kg s' 1 . Total NO x mass is given in 10 5 kg. 

Figure 5: As Figure 3, but for ozone (ppbv). 

Figure 6: Frequency distributions for (left) and ozone (right) derived from NOXAR 
measurements (top), ECHAM4.L39(DLR)/CHEM data (middle) and GISS data (bottom). 
Model data represent 12 hourly means. 

Figure 7: Contributions [in % of the control simulation] of aircraft (red), lightning (green) and 
surface (yellow) NO x emissions as well as the upper boundary (blue) upper troposphere and 
lower stratosphere NO x (right) and ozone (left) concentration derived with the 
ECHAM4.L39(DLR)/CHEM (top) and GISS-model (bottom). Values are given for DJF 
(solid) and JJA (dotted). For the GISS model values above the tropopause are defined by the 
upper boundary and therefore either 0% or 100%, indicated with thinner lines. 

Figure 8: Frequency distributions for NO x (top) and ozone (bottom) derived with the 
ECHAM4.L39(DLR)/CHEM and GISS-model for the control runs (solid) and perturbation 
runs without lightning (long dashed), aircraft (short dashed), surface emissions (dotted) and 
without a stratospheric flux (dotted-dashed). 
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